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Satellite  Observations  of  Marine  Stratus/Stratocumulus 
ABSTRACT 

Work  supported  by  this  ONR  Grant  led  to  the  stream  lining  of  spatial  coherence  analysis 
software  and  allowed  the  application  of  the  spatial  coherence  analysis  to  a  large  volume  of 
satellite  imagery  data.  The  analyses  were  performed  in  order  to  obtain  the  properties  of  low-level 
clouds  and  the  cloud-ffee  ocean  background.  The  work  led  to  the  findings  that;  1)  reflectivities  of 
broken  clouds  at  both  0.63  and  3.7  pm  are  smaller  than  those  of  nearby  untform  clouds  indicating 
that  the  broken  clouds  have  substantially  less  liquid  water  than  do  their  uniform  counteqjarts  and 
that  photons  are  escaping  through  the  sides  of  the  broken  clouds  and  being  absorbed  by  the  ocean 
surface.  2)  The  thin  atmosphere  approximation  used  to  justify  the  split  window  technique  for 
obtaining  sea  surface  temperatures  from  satellite  imagery  data  fails  to  give  the  sensitivity  of  the 
emitted  radiances  to  variations  in  the  atmospheric  optical  path  that  arise  from  changes  in  the  scan 
angle  of  the  imaging  radiometer.  3)  Substantial  portions  of  low-level  marine  stratocumulus  are 
semitransparent  at  1 1-pm.  4)  As  do  fractional  cloud  cover  and  liquid  water  amount,  droplet 
radius  appears  to  vary  from  pixel  to  pixel  at  the  scale  available  from  conventional  satellite  imagery 
data  (1  km).  5)  Comparisons  of  cloud  cover  derived  using  the  spatial  coherence  method  applied 
to  satellite  imagery  data  for  the  Atlantic  Stratocumulus  Transition  Experiment  (ASTEX)  with 
cloud  cover  predicted  by  forecast  models  indicated  that  the  forecast  models  were  unable  to 
adequately  predict  low-level  clouds.  The  deficiency  was  thought  to  be  due  to  the  inability  of  the 
model  to  adequately  predict  the  relative  humidities  in  the  upper  levels  of  the  boundary  layer. 
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1. 


Introduction 


The  goal  of  this  project  was  to  develop  automated  procedures  incorporating  the  spatial 
coherence  method  for  analyzing  cloud  properties  from  satellite  imagery  data  and  then  to  use  the 
analysis  scheme  to  analyze  the  cloud  properties  obtained  during  the  Atlantic  Stratocumulus 
Transition  Experiment  (ASTEX).  ASTEX  took  place  in  the  Northeastern  Atlantic  in  June  1992. 
During  the  course  of  the  project,  the  spatial  coherence  analysis  software  was  streamlined  and 
large  volumes  of  satellite  imagery  data  were  analyzed.  The  analyzed  results  afforded 
opportunities  to  characterize  the  properties  of  marine  stratocumulus  and  to  test  the  assumptions 
underlying  the  estimation  of  sea  surface  temperatures  (SSTs)  from  multispectral  infrared  imagery 
data.  Also,  a  new  approach  was  developed  for  analyzing  cloud  properties  from  multispectral 
infrared  imagery  data.  The  new  approach  provided  a  test  of  the  opaque  cloud  approximation  used 
in  the  spatial  coherence  method  to  obtain  the  fractional  cloud  cover  for  imager  pixels  (~1  —  4  km) 
that  are  only  partly  cloud  covered.  The  new  approach  itself  was  subjected  to  consistency  checks 
with  the  1-km  Advanced  Very  High  Resolution  Radiometer  (AVHRR)  data  collected  during 
ASTEX.  Finally,  the  cloud  cover  derived  from  the  ASTEX  observations  using  the  spatial 
coherence  method  was  used  to  assess  the  performance  of  cloud  parameterization  schemes  used  in 
weather  forecast  models.  These  various  projects  and  their  findings  are  described  in  the  following 
sections. 

2.  Reflectivities  of  Uniform  and  Broken  Layered  Cloud  Systems 

The  spatial  coherence  method  was  used  to  obtain  cloud  properties  for  ASTEX.  It  was 
applied  to  the  1-km  AVHRR  data  collected  by  the  Naval  Post  Graduate  School  during  the 
experiment.  The  special  feature  of  the  spatial  coherence  method  is  that  it  distinguishes  between 
imager  pixels  that  are  cloud-free,  overcast,  and  partly  cloud-cloud  covered.  Following  an  earlier 
study  (Coakley  ,  1991),  the  data  collected  during  ASTEX  was  analyzed  to  obtain  the  reflectivities 
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of  low-level,  layered  cloud  systems  at  0.63  and  3.7  pm.  Reflectivities  of  low-level  clouds  that 
completely  covered  several  of  the  1-km  pixels  were  obtained  along  with  the  reflectivities  of  clouds 
that  only  partly  covered  the  1-km  pixel.  The  spatial  coherence  results  were  used  to  segregate 
overcast  pixels  from  partly  cloudy  pixels.  Comparisons  were  made  between  the  reflectivities  for 
the  overcast  and  broken  clouds  within  localized  (~  60  km)  regions  so  that  the  properties  of  the 
uniform  and  broken  clouds  were  arguably  identical  save  for  their  different  extents  cloud  cover. 

Figure  1  shows  the  reflectivities  of  low-level,  layered  clouds  at  0.63  and  3.7  pm  derived 
from  the  overcast  pixels  (solid  curves)  and  differences  in  the  reflectivities  (overcast  —  broken) 
(dashed  curves)  obtained  for  the  nearby  partly  cloudy  pixels  collected  during  ASTEX.  The 
observations  are  for  the  afternoon  passes  of  the  NOAA-11  satellite.  The  reflectivities  are 
presented  as  functions  of  satellite  zenith  angle.  Positive  zenith  angles  indicate  radiation  that  is 
reflected  in  the  direction  of  forward  scattering;  negative  zenith  angles  indicate  radiation  that  is 
reflected  in  the  direction  of  backward  scattering.  Points  indicate  mean  values  and  error  bars 
indicate  the  expected  uncertainty  (one  standard  deviation)  in  the  mean  values.  The  uncertainty  is 
estimated  on  the  basis  of  the  variability  of  the  observations  obtained  during  ASTEX. 

Like  the  results  reported  by  Coakley  and  Judge  (1991)  for  a  set  of  pre-ASTEX 
observations,  the  ASTEX  observations  conform  to  the  findings  reported  by  Coakley  (1991).  At 
the  visible  wavelength,  broken  clouds  have  lower  reflectivities  than  their  uniform  counterparts. 

As  with  the  earlier  study,  the  discount  factor  is  approximately  0.8.  This  reduction  is  due  either  to 
smaller  column  amounts  of  cloud  liquid  water  in  the  broken  clouds  or  to  the  escape  of  photons 
from  the  sides  of  the  broken  clouds  and  their  subsequent  absorption  by  the  ocean  surface. 
Regardless  of  the  mechanism,  the  anisotropy  of  the  radiation  reflected  by  broken  clouds  is,  within 
the  uncertainty  of  the  observations,  similar  to  that  reflected  by  uniform  clouds. 
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At  3.7  pm,  reflectivities  of  broken  clouds  are  smaller  than  those  of  their  uniform 
counterparts  for  radiation  that  is  reflected  in  the  direction  of  backward  scattering  (indicated  by 
negative  satellite  zenith  angles)  and  larger  than  the  reflectivities  of  their  uniform  counterparts  for 
radiation  that  is  reflected  in  the  direction  of  forward  scattering  (indicated  by  positive  satellite 
zenith  angles).  The  higher  reflectivities  for  broken  clouds  in  the  forward  scattering  direction 
should,  however,  be  viewed  with  caution.  At  3.7  pm,  cloud-free  oceans  exhibited  large 
reflectivities  (not  shown),  even  larger  than  those  for  overcast  pixels.  The  high  reflectivities 
correspond  with  high  reflectivities  at  0.63  pm  and  are  evidently  due  to  sunglint.  Because  of  the 
sunglint,  estimates  of  reflectivities  for  clouds  that  only  partly  cover  the  pixels  may  be  subject  to 
large  errors.  Differences  in  the  reflectivities  for  the  radiation  reflected  in  the  direction  of 
backward  scattering,  however,  are  reliable.  As  with  the  reflectivities  at  visible  wavelengths,  the 
reduction  in  reflectivities  for  the  broken  clouds  can  be  explained  either  due  to  a  decrease  in  the 
amount  of  cloud  water  for  broken  clouds  or  to  the  escape  of  photons  from  the  sides  of  the  clouds. 
Unlike  the  reflectivities  at  visible  wavelengths,  however,  the  anisotropy  of  3.7-pm  radiation 
reflected  by  broken  clouds  would  appear  to  differ  from  that  reflected  by  uniform  clouds. 

Radiation  reflected  by  broken  clouds  would  be  more  isotropic  than  that  reflected  by  uniform 
clouds. 


3.  SSTs  from  AVHRR  Observations 

The  spatial  coherence  analysis  was  applied  to  4-km  AVHRR  observations  for  the  months 
of  May,  June,  and  July  1983  and  1984  to  determine  the  types  of  cloud  conditions  that  were  likely 
to  be  found  during  ASTEX.  In  addition  to  the  information  on  cloud  properties,  these 
observations  provided  opportunities  to  test  concepts  commonly  used  for  retrieving  sea  surface 
temperatures  (SSTs)  from  1 1  and  12-pm  imagery  data.  As  part  of  the  work  undertaken  under 


4 


this  ONR  grant,  Mr.  A.  Kowalski  completed  the  requirements  for  an  MS  degree  in  Atmospheric 
Sciences  at  OSU. 

Kowalski  explored  the  possibility  of  inferring  SSTs  from  the  relationships  exhibited  by  the 
emitted  radiances  themselves,  i.e.  without  the  need  to  obtain  collocated  buoy  or  shipboard 
measurements  of  SSTs  (Kowalski,  1993).  The  approach  was  based  on  the  analysis  of  the  split- 
window  technique  for  obtaining  SSTs  given  by  McMUlin  (1975)  and  by  McMilhn  and  Crosby 
(1984).  According  to  McMillin,  the  radiance  observed  in  channel  4  (11  pm)  of  the  AVHRR 
should  be  linearly  related  to  the  difference  in  the  radiances  observed  in  channels  4  and  5  (12  pm). 
The  relationship  is  given  by 


^4  ~  ^4(^5)  ■*'7  (^4  ^5) 

where  £4(7^)  is  the  Planck  function  at  the  surface  temperature  ,  which  is  the  desired  SST,  y  is 

a  constant,  and  is  the  radiance  in  channel  4  that  would  be  observed  if  the  observed  channel  5 

brightness  temperature  were  used  to  calculate  the  channel  4  radiance.  Usually,  y  is  determined 

by  correlating  (1)  with  in  situ  buoy  or  shipboard  observations  of  .  The  form  of  (1),  however, 

suggests  that  if  a  region  has  a  sufficiently  uniform  surface  temperature  so  that  is  constant, 

/ 

then  y  is  the  slope  for  the  relationship  between  and 

Kowalski  searched  5°  x  5°  latitude  -  longitude  regions  containing  numerous  cloud-free 

pixels  as  determined  by  spatial  coherence  analysis  of  4-km  AVHRR  observations  for  May- July 

1983  and  1984  in  the  north  Adantic  (5°N  -  45°N,  10°W  -  40°W).  Few  regions  were  found  in 

/ 

which  a  linear  relationship  between  and  -  /j  could  be  claimed.  The  lack  of  linear 
relationships  was  probably  caused  by  sea  surface  temperature  gradients  within  the  ~  500  km 

regions.  The  analysis  of  smaller  regions  was  impractical  because  significant  variations  of  and 
/ 

-/j  across  the  region  are  needed  in  order  to  obtain  a  reliable  estimate  of  y  .  If  the  SST  is 
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uniform  across  the  region,  variations  of  and  — /j  are  due  to  changes  in  the  atmospheric 

optical  path  which  come  with  changes  in  satellite  zenith  angle  (i.e.  with  changes  in  the  scan  angle 
oftheAVHRR). 


Figure  2  shows  one  incident  in  which  a  relationship  in  accordance  with  (1)  was  observed. 
In  this  case  y  =  -  4.8.  The  commonly  used  value  is  y  ~  -  3.  The  magnitude  of  y  increases  as 
the  burden  of  water  vapor  in  the  atmosphere  increases. 

The  lack  of  regions  with  sufficient  uniformity  in  SST  to  obtain  values  of  y  was  surprising 
and  led  to  the  assessment  of  the  “thin  atmosphere”  approximation  that  gives  rise  to  (1).  Using  a 
radiative  transfer  model  for  the  emitted  radiances  observed  by  the  AVHRR,  Kowalski  calculated 
radiances  as  a  function  of  satellite  zenith  angle  for  tropical,  midlatitude  and  polar  profiles  of 
temperature  and  moisture.  The  results  are  shown  in  Figure  3.  The  figure  demonstrates  that  (1)  is 
unsuitable  for  changes  in  the  radiance  due  to  the  changes  in  the  optical  path  that  come  with 
changes  in  satellite  zenith  angle.  Owing  to  the  nonlinear  relationship  between  the  emitted 
radiances  and  the  column  amount  of  water  vapor,  the  slopes  obtained  on  the  basis  of  changes  in 
the  radiances  with  satellite  zenith  angle  (dashed  lines)  are  grossly  biased  values  of  y  (solid  lines) 
for  tropical  and  midlatitude  profiles.  The  slope  is  a  useful  estimate  only  in  the  case  of  a  polar 
profile,  which  has  little  water  vapor,  and  consequently,  is  appropriate  for  the  “thin  atmosphere” 
approximation. 


4.  Retrieval  of  Cloud  Properties  from  Multi-Spectral  Infrared  Imagery  Data 

Clouds,  particularly  low-level  maritime  clouds,  cannot  be  expected  to  fill  the  field  of  view 
of  a  radiometer  having  even  a  1-km  resolution.  Broken  clouds  within  a  field  of  view  should  be 
the  norm.  Clouds,  on  the  other  hand,  usually  come  in  well-defined  layers.  By  finding 
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geographically  localized  regions  (~  4  —  8  km)  of  uniform  emission,  the  spatial  coherence  method 
utilizes  the  tendency  of  clouds  to  form  layers  to  identify  pixels  which  are  overcast  by  layered 
clouds.  The  localized  geographic  regions  which  lack  uniform  emission  are  presumed  to  contain 
broken  clouds.  In  order  to  obtain  the  fractional  cloud  cover  of  regions  with  broken  clouds,  the 
assumption  is  made  that  the  clouds  within  the  region  being  analyzed  are  opaque  and  part  of  the 
layer,  or  layers,  indicated  by  the  localized  regions  exhibiting  uniform  emission.  Through  NASA 
support  for  the  FIRE  project,  Lin  and  Coakley  (1993)  developed  a  new  approach  for  obtaining 
the  properties  of  layered  cloud  systems  from  imagery  at  two  infrared  wavelengths.  Work 
supported  in  part  through  this  ONR  grant  used  the  Lin  and  Coakley  retrieval  scheme  to 
characterize  the  properties  of  marine  stratocumulus  off  the  coast  of  Chile  (Luo  et  al.,  1994). 
Work  also  supported  in  part  through  this  ONR  grant  extended  the  Lin  and  Coakley  retrieval 
scheme  to  determine  whether  there  was  consistency  in  the  retrievals  when  all  five  channels  of  the 
AVHRR  were  utilized  (Coakley  and  Lin,  1993). 

For  single-layered  cloud  systems,  the  emitted  radiance  observed  by  channel  i  of  a  satellite 
radiometer  is  given  by 


A  =  (1  ~  A)  hi  +  Ai^ihi + 


(2) 


where  is  the  radiance  associated  with  the  cloud-free  portion  of  the  pixel,  e,  is  the  average 
emissivity  of  the  clouds  in  the  pixel  at  the  wavelengths  associated  with  channel  /,  is  the 

emission  that  would  be  observed  if  the  pixels  were  overcast  by  the  clouds  and  if  the  clouds  were 
black  at  the  wavelengths  associated  with  channel  i  and  t-  is  the  average  transmissivity  of  the 

clouds.  Here  it  is  assumed  that  the  observations  are  being  made  in  an  atmospheric  "window"  so 
that  emission  by  the  atmospheric  gases  above  the  clouds  is  taken  to  be  negligible.  As  is  illustrated 
in  Figure  4,  based  on  (2),  emitted  radiances  from  two  distinct  spectral  intervals,  (e.g.,  11  and 
12  pm  or  3.7  and  12  pm)  are  constrained  to  fill  the  radiance  domain  bracketed  by  well-prescribed 
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limits.  If  the  clouds  are  opaque,  then  =  0  and  e,  -  where  6,max  “  ^  ^max  ^max  the 
maximum  value  of  the  reflectivity  that  can  be  obtained  at  the  wavelengths  associated  with  channel 
i.  At  some  wavelengths  where  absorption  is  appreciable,  like  at  1 1  and  12  p,m,  taken 

to  be  zero.  According  to  (2)  when  the  clouds  are  part  of  a  single  layer  and  they  are  opaque  at  any 
two  thermal  infrared  wavelengths,  then  the  radiances  at  these  wavelengths  are  linearly  related. 
Such  relationships  have  often  been  observed  at  night  for  3.7  and  1 1-lim  radiances  (Coakley  and 
Davies,  1986).  For  pixels  that  are  overcast  by  single-layered  clouds  that  are  semitransparent  at 
the  infrared  wavelengths  associated  with  one  or  both  channels,  then  again  according  to  (2),  the 
radiances  will  be  nonlinearly  related  through  the  separate  nonlinear  dependencies  of  emissivity  and 
transmissivity  on  the  optical  path  length.  Finally,  for  pixels  that  are  overcast  by  clouds  that  are 
opaque  for  any  two  infrared  channels  and  are  distributed  in  altitude,  the  radiances  will  follow, 
aside  from  effects  due  to  absorption  and  emission  in  the  cloud-free  atmosphere,  the  Planck 
blackbody  curve. 


Using  various  values  of  in  (2)  along  with  various  values  for  the  liquid/ice  water  column 
amounts  for  droplets  with  a  fixed  size  distribution  to  calculate  values  of  e,  and  r, ,  a  set  of  curves 
are  generated  which  may  be  used  to  interpret  the  observed  radiances.  The  curves  are  shown  in 
Figure  4.  For  60  km  regions,  single-layered  cloud  systems  are  often  broken.  They  often  have 
radiative  properties  which  are  consistent  with  a  range  of  optical  pathlengths.  Consequently, 
instead  of  giving  rise  to  radiances  that  fall  along  a  particular  line,  whether  it  be  linear  or  nonlinear, 
the  radiances  appear  to  fill  the  radiance  domain  bracketed  by  the  extremes.  An  example  is  shown 
in  Figure  5.  Cloud  cover  and  optical  pathlength  (or  equivalently,  as  is  given  in  the  figure,  the 
1  l-|im  emissivity)  are  then  assigned  to  a  pixel-scale  radiance  pair  on  the  basis  of  the  radiances 
predicted  by  the  radiative  transfer  model.  The  results  shown  in  the  figure  are  obtained  using  the 
Eddington  approximation  to  calculate  e,  and  . 
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As  was  noted  by  Lin  and  Coakley  (1993),  the  size  of  the  domain  covered  by  the  infrared 
radiances  is  a  function  of  hydrometeor  size.  The  area  covered  is  generally  large  for  small  droplets 
and  collapses  to  a  straight  line  (no  area)  for  large  droplets.  Cloud  properties  are  retrieved  from 
radiance  pairs  by  adjusting  the  effective  particle  radius  (and  also  the  emission  associated  with 
opaque,  overcast  clouds)  in  the  radiative  transfer  model  until  the  predicted  radiance  domain 
tightly  hugs  the  observed  radiances  (Lin  and  Coakley,  1993).  An  example  is  shown  in  Fig.  5. 

Luo  et  al.  (1993)  used  4-km  resolution  emitted  radiances  at  11  and  12  |j.m  from  the 
NOAA-9  AVHRR  to  survey  the  properties  of  marine  stratocumulus  off  the  coast  of  South 
America.  They  found  droplet  radii  that  vary  within  the  typically  observed  ranges  (-8-12  pm). 
One  of  the  surprising  results  of  this  study  was  that  a  substantial  fraction  of  the  area  covered  by 
marine  stratocumulus  is  semitransparent  at  11  pm.  The  250  km  average  11 -pm  emissivity  was 
between  0.7  and  0.8.  These  results  were  consistent  with  the  earlier  findings  of  Wielicki  and 
Parker  (1992).  Conventional  wisdom  took  marine  stratocumulus  to  be  black.  Radiative  transfer 
and  aircraft  observations  had  indicated  that  such  clouds  become  opaque  when  they  reach  depths 
of  100  -  450  m,  a  common  occurrence  (Platt,  1976).  But  during  field  experiments,  thin  clouds  are 
often  avoided  as  the  mission  is  generally  to  obtain  the  properties  of  clouds:  so,  the  thicker  the 
cloud  the  better.  As  a  result,  there  is  abundant  evidence  from  field  programs  that  such  clouds  are 
black. 


Not  only  were  marine  stratocumulus  found  to  be  semitransparent,  but  also  their  average 
emissivity  was  found  to  be  a  function  of  the  regional  scale  cloud  cover  as  is  shown  in  Figure  6. 
Evidently,  if  the  region  is  heavily  cloud  covered,  the  ratio  of  semitransparent  to  opaque  cloud 
material  decreases.  The  consequence  is  that,  not  surprisingly,  regional  scale  cloud  cover  and 
liquid  water  column  amounts  are  related  to  each  other. 
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If  this  multispectral  infrared  retrieval  procedure  is  to  be  believed,  then  cloud  cover, 
liquid/ice  water  column  amounts  and  particle  size  obtained  using  one  set  of  wavelengths,  e.g.  11 
and  12  |J,m  should  match  those  obtained  using  a  second  set,  e.g.  3.7  and  12  [im  (Coakley  and  Lin 
1993).  There  should  be  a  one-to-one  mapping  of  the  pixels  in  the  two-channel  radiance  domains 
as  is  illustrated  in  Figure  7.  Figures  8  and  9  show  an  experiment  that  tests  this  mapping. 

Nighttime  observations  are  used  in  order  to  avoid  the  reflection  of  sunlight  at  3.7  fim.  In  the  test 
the  11  and  12-|im  radiances  are  used  to  identify  clusters  of  pixels  that  would  appear  to  be 
overcast  by  semitransparent  clouds  (squares)  and  clusters  of  pixels  that  would  appear  to  be  partly 
covered  by  opaque  clouds  (triangles)  (Fig.  8).  The  location  of  the  pixel  clusters  in  the  3.7-12-|xm 
radiance  domain  is  shown  in  Fig.  9.  Clearly,  the  interpretation  based  on  the  1 1-12-p.m  radiance 
pairs  is  inconsistent  with  that  based  on  the  3.7-12-|im  radiance  pairs. 

Even  though  only  one  example  is  presented,  the  example  is  representative  of  observations 
for  single-layered  cloud  systems  found  at  night.  The  same  principles  were  used  to  test  the 
consistency  of  emitted  radiation  at  1 1  and  12  pm  with  emitted  radiation  at  1 1  pm  and  reflected 
sunlight  at  0.63  pm.  As  with  emission  at  3.7, 11  and  12  ^im,  the  interpretation  based  on  emission 
at  1 1  and  12  |j.m  does  not  hold  for  observations  at  0.63  and  1 1  iim.  The  conclusion  drawn  from 
the  results  of  these  tests  is  that  effective  particle  size,  like  fractional  cloud  cover  and  cloud  optical 
path,  varies  from  pixel  to  pixel  at  the  1-km  spatial  resolution  available  to  the  AVHRR.  A  simple 
solution  to  rectify  this  inconsistency  would  be  to  utilize  observations  at  say,  3.7, 11  and  12  pm  to 
derive  fractional  cloud  cover,  particle  size  and  cloud  optical  path  on  a  pixel  by  pixel  basis.  Arking 
and  Childs  (1985)  proposed  such  a  scheme,  albeit  using  reflected  sunlight  at  0.63  pm  and  3.7  p,m 
and  emitted  radiances  at  1 1  pm.  Work  on  this  problem  is  continuing  through  support  by  NASA. 


10 


5. 


Evaluation  of  ECMWF  Cloud  Parameterization  Scheme  during  ASTEX 


Spatial  coherence  analyses  of  cloud  cover  and  cloud  top  height  (as  indicated  by  the 
observed  1 1-pm  brightness  temperature  for  overcast  pixels)  during  ASTEX  were  compared  with 
surface  based  estimates  of  cloud  cover  from  cielometer  data  and  predictions  of  cloud  cover  using 
ECMWF  model  forecast  and  analysis  fields  combined  with  various  diagnostic  cloud 
parameterization  schemes  (Bretherton  et  al.,  1995).  Cloud  properties  obtained  with  the  spatial 
coherence  method  applied  to  1-km  AVHRR  data  were  analyzed  on  a  2  x  2  latitude  —  longitude 
grid.  At  this  scale,  low-level  cloud  systems  were  often  accompanied  by  upper  level  cloud 
systems,  impeding  the  analysis  of  conditions  when  low-level  clouds  existed  by  themselves.  Figure 
10  shows  a  comparison  of  the  analyzed  satellite  derived  fractional  cloud  cover  with  that  inferred 
for  the  corresponding  three-hour  period  of  cielometer  observations  on  Porto  Santo  when  only 
low-level  clouds,  as  deduced  from  the  satellite  observations  were  present.  The  agreement  is 
satisfying,  if  somewhat  surprising,  given  the  differences  in  the  spatial  scales  associated  with  the 
two  different  observations.  The  cielometer  observations  are  for  clouds  directly  over  the  site. 

With  the  7  m  s~‘  surface  wind  observed  during  ASTEX,  the  average  for  the  three  hour  period 
would  seem  appropriate  for  a  ~  70  km  region  downwind  of  the  site.  The  satellite  analysis,  on  the 
other  hand,  is  for  200  km  scale  regions.  Nevertheless,  high  values  of  the  cloud  cover  in  the 
satellite  observations  correspond  with  high  values  in  the  cielometer  observations  and  conversely 
low  values  in  the  satellite  observations  correspond  with  low  values  in  the  cielometer  observations. 


The  agreement  between  satellite  and  cielometer  based  observations  is  to  be  contrasted 
with  the  lack  of  agreement  between  the  cloud  cover  obtained  with  the  ECMWF  cloud 
parameterization  scheme  and  the  cielometer  cloud  cover  shown  in  Figure  1 1.  The  cloud  cover 
predicted  using  the  ECMWF  scheme  with  either  the  model  generated  or  analyzed  fields  was  too 
low.  The  chief  deficiency  in  the  model  was  thought  to  be  the  relatively  low  values  of  the  relative 
humidity  near  the  top  of  the  boundary  layer  produced  by  the  model  (Bretherton  et  al.,  1995). 
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6. 


Summary  of  Findings 


The  work  supported  under  this  ONR  grant  led  to  the  following  findings: 

•  Reflectivities  at  0.63  pm  and  3.7  pm  for  broken,  low-level,  layered  clouds  observed 
during  ASTEX  were  smaller  than  the  reflectivities  of  their  uniform  counterparts.  The 
lower  values  may  be  explained  by  either  lower  amounts  of  cloud  liquid  water  in  broken 
clouds  or  by  the  escape  of  photons  through  the  sides  of  broken  clouds. 

•  The  optically  thin  approximation  used  to  justify  split  window  techniques  for  the  estimation 
of  sea  surface  temperature  from  satellite  infrared  imagery  data  cannot  be  used  to  infer  the 
variation  of  the  emitted  radiances  with  satellite  zenith  angle.  The  relationship  between  the 
emitted  radiances  and  the  atmospheric  optical  pathlength  is  nonlinear  and  outside  of  the 
range  appropriate  for  the  optically  thin  approximation  to  apply. 

•  A  substantial  portion  of  marine  stratocumulus  resides  in  cloud  material  that  is 
semitransparent  at  11  pm.  The  average  11 -pm  cloud  emissivity  for  250-km  scale  regions 
lies  between  0.7  and  0.8. 

•  Like  fractional  cloud  cover  and  liquid  water  path,  droplet  radius  for  marine  stratus  appears 
to  vary  from  pixel  to  pixel  at  the  highest  resolutions  available  to  conventional  satellite 
observations  (~  1  km). 
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Owing  to  their  inability  to  predict  the  high  levels  of  relative  humidity  which  appear  to  be 
present  in  the  upper  levels  of  the  boundary  layer,  forecast  models  are  unable  to  predict  the 
presence  of  low-level  cloud. 
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Figure  Captions 


Figure  1.  Reflectivities  for  overcast  pixels  (solid  line)  and  differences  in  reflectivities  (overcast  - 
broken)  at  0.63  pm  (a)  and  3.7  pm  (b)  obtained  for  low-level,  layered  cloud  systems 
observed  during  ASTEX.  The  symbols  indicate  average  values  for  the  indicated  satellite 
zenith  angle;  the  error  bars  indicate  estimates  of  the  uncertainty  (one  standard  deviation) 
in  the  mean  values.  Positive  satellite  zenith  angles  indicate  the  forward  scattering 
direction;  negative  values  indicate  the  backward  scattering  direction. 

Figure  2.  Linear  relationship  between  the  observed  cloud-free  radiance  in  channel  4(11  pm), 
of  the  AVHRR  and  the  difference  between  the  channel  4  radiance  and  the  channel  4 

/ 

radiance  that  would  be  obtained  for  the  channel  5  (12  pm)  brightness  temperature 
The  observations  are  for  a  5°  x  5°  latitude  -  longitude  region.  The  linear  relationship  is 
expected  if  the  sea  surface  temperature  is  uniform  across  the  region.  The  variation  of 
and  /j  across  the  region  is  presumed  to  be  due  to  the  variation  of  optical  path  with 
AVHRR  scan  angle  across  the  region. 

Figure  3.  Slopes  giving  y  in  (1).  Solid  lines  are  the  actual  values  of  y  .  Dashed  lines  are  the 

/ 

values  of  y  predicted  on  the  basis  of  the  variations  in  and  with  satellite  zenith 

angle. 

Figure  4.  11  and  12-pm  radiances  given  by  (2)  for  a  single-layered,  low-level  cloud  system 

exhibiting  various  degrees  of  fractional  cloud  cover  and  cloud  water  path  as  indicated  by 
the  11-pm  emissivities.  Solid  lines  are  lines  of  constant  cloud  cover  with  the  cloud  cover 
given  by  the  numbers  above  the  lines.  Dashed  lines  are  lines  of  constant  1 1-pm  emissivity 
with  the  emissivities  given  at  the  ends  of  the  lines.  The  calculations  were  performed  for  a 
cloud  with  mean  droplet  radius,  R^fi  =  2.7  pm. 
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Figure  5.  Observed  and  theoretically  calculated  1 1  and  12-pm  radiances. 


Figure  6.  Average  cloud  cover  and  1 1-pm  emissivity  for  250-km  regions  containing  low-level 
marine  stratocumulus  off  the  coast  of  Chile. 

Figure  7.  Bispectral  consistency  which  should  exist  if  the  assumptions  underlying  the 
multispectral  retrieval  procedures  are  valid. 

Figure  8.  Identification  of  pixel  clusters  containing  opaque,  broken  clouds  (triangles)  and 
overcast  semitransparent  clouds  (squares)  based  on  emission  at  1 1  and  12  pm. 

Figure  9.  Mapping  of  pixel  clusters  containing  opaque,  broken  clouds  (triangles)  and  overcast 

semitransparent  cloud  (squares)  as  identified  through  their  emission  at  11  and  12  pm.  For 
a  consistent  mapping,  the  squares  should  fall  on  the  upper  envelope  of  the  3.7- 12-pm 
radiance  domain  occupied  by  the  radiance  pairs  and  the  triangles  should  fall  along  the 
straight  line. 

Figure  10.  Spatial  coherence  estimates  of  low-level  cloud  cover  and  corresponding  estimates 
based  on  ceilometer  observations  at  Porto  Santo  during  the  ASTEX  experiment. 

Figure  11.  ECMWF  estimates  of  low  level  cloud  cover  and  corresponding  estimates  based  on 
ceilometer  observations  at  Porto  Santo  during  the  ASTEX  experiment. 
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11-12  Radiance  Pairs 


LOCAL  MEAN  1  1  /Li.m  RADIANCE 


1 .  Single-Layer,  plane-parallel. 

2.  Eddington  approximation. 

3.  Observed  cloud-free  and  overcast  radiances. 

4.  Homogeneous  droplet  radii  given  by  log-normal 
size  distribution. 


Figure  4 


Fit  to  Observations 
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2x2  arrays  of  1-km  HRPT  data  from  nighttime  overpass  of  NOAA-11 
for  60-km  scale  portion  of  ASTEX  region  containing  single-layered  marine 
stratocumulus  system.  10  June  1992,  0412  Z 
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Figure  6 
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Bispectral  Consistency 
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Figure  7 
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Pixel  Cluster  Locations  in  the 
Multispectral  Radiance  Domain 
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Triangles— pixels  containing  opaque,  broken  cloud. 
Squares— pixels  overcast  by  semitransparent  cloud 

2x2  arrays  of  1-km  HRPT  data  from  nighttime  overpass  of 
NOAA-11  for  60-km  scale  portion  of  ASTEX  region  containing 
single-layered  marine  stratocumulus  system.  10  June  1992, 
0412  Z 


Figure  8 
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Pixel  Cluster  Locations  in  the 
Multispectral  Radiance  Domain 

12-3.7  fivn  Radiance  Pairs 


Triangles— pixels  containing  opaque,  broken  cloud. 
Squares— pixels  overcast  by  semitransparent  cloud 

2x2  arrays  of  1-km  HRPT  data  from  nighttime  overpass  of 
NOAA-11  for  60-km  scale  portion  of  ASJEX  region  containing 
single-layered  marine  stratocumulus  system.  10  June  1992, 
0412  Z 


Figure  9 
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Porto  Santo  low  cloud: 

ECMWF  vs.  3  hr  avg  ceilometer  cloud  amount  correlation 
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